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We predict (7^!°^^, the proton-air inelastic production cross section, at pp center-of-mass energies 
2 < y/s < 100000 TeV, using high energy predictions from a saturated Froissart bound parame- 
terization of accelerator data on forward pp and pp scattering amplitudes, together with Glauber 
theory. The parameterization of the pp and pp cross sections incorporates analyticity constraints and 
unitarity, allowing accurate extrapolations to ultra-high energies. Our predictions are in excellent 
agreement with cosmic ray extensive air shower measurements, both in magnitude and in energy 
dependence. 

PACS numbers: 13.60.Hb, 12.38.-t, 12.38.Qk 



Introduction. There are now available published p- 
air inelastic production cross sections [J, 0, H, 0] (crp™^^.) 
that span the enormous pp cms (center-of-mass system) 
energy range 2 < ^/s < 100000 TeV, reaching energies 
well above the Large Hadron Collider (LHC). Moreover, 
there are also now available very accurate predictions at 
cosmic ray energies for the total pp cross section, app, 
from fits[5| to accelerator data that used adaptive data 
sifting algorithms fi*! and analyticity constraints^. How- 
ever, extracting proton-proton cross sections from pub- 
lished cosmic ray observations of extensive air showers, 
and vice versa, is far from straightforward 8] . By a vari- 
ety of experimental techniques, cosmic ray experiments 
map the atmospheric depth at which extensive air show- 
ers develop and measure the distribution of Xmaxj the 
shower maximum, which is sensitive to the inelastic p-air 
cross section afj°!^:^. From the measured X,nax distribu- 



p~air 



tion, the experimenters deduce fp^^r- ^^^^ note we 
will compare published values of Up^^j. with predictions 



made from app, using a Glauber model to obtain a, 



pp- 



from a 

prod 
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'''p-air from the Xtnax distribution: Method I. The mea- 



sured shower attenuation length (A^) is not only sensi- 
tive to the interaction length of the protons in the atmo- 
sphere (Ap_air), with 



A™ — kXr, 



prod 
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= k 



24, 100 

prod 



(1) 



(with Am and Ap_air in gem ^, the proton mass m in g, 
and the inelastic production cross section crpL°air in mb), 
but also depends on the rate at which the energy of the 
primary proton is dissipated into electromagnetic shower 
energy observed in the experiment. The latter effect is 
parameterized in Eq. ([1]) by the parameter k. The value 
of k depends critically on the inclusive particle produc- 
tion cross section and its energy dependence in nucleon 
and meson interactions on the light nuclear target of the 
atmosphere (see Ref . Q ) . We emphasize that the goal of 
the cosmic ray experiments is a-p!°air (^^ correspondingly, 
Ap-air), whereas in Method I, the measured quantity is 



TABLE I: A table of fc-values, used in experiments and from 
Monte Carlo model simulation 



Experiment 


k 


Fly's Eye 


1.6 


AGASSA 


1.5 


Yakutsk 


1.4 


EASTOP 


1.15 


Monte Garlo Results: C.L. Pryke 


Model 


k 


GORSIKA-SIBYLL 


1.15 ±0.05 


MOCCA-SIBYLL 


1.16 ±0.03 


GORSIKA-QGSjet 


1.30 ±0.04 


MOCCA-Internal 


1.32 ±0.03 



A„i- Thus, a significant drawback of Method I is that 
one needs a model of proton-air interactions to complete 
the loop between the measured attenuation length Am 
and the cross section CTp^firJ ^-^v one needs the value of 
k in Eq. ^ to compute cr^^^j. . Shown in Table [J are 
the widely varying values of k used in the different ex- 
periments. Clearly the large range of fc-values, from 1.15 
for EASTOPQ to 1.6 for Fly's EyeQ differ significantly, 
thus making the published values of crpl"^^ unreliable. It 
is interesting to note the monotonic decrease over time in 
the fc's used in the different experiments, from 1.6 used 
in Fly's Eye in 1984 to the 1.15 value used in EASTOP in 
2007, showing the time evolution of Monte Carlo models 
of energy dissipation in showers. For comparison, Monte 
Carlo simulations made by PrykeQ in 2001 of several 
more modern shower models are also shown in Table 
m We see that even among modern shower models, the 
spread is still significant. The purpose of this letter is a 
proposal to minimize the impact of model dependence on 
the determination of o'p'_°J^j. 

(Tp™^j. from the ATmax distribution: Method II. The 
HiRes group [l3| has developed a quasi model- free method 
of measuring cTp^^j. directly. They fold into their shower 
development program a randomly generated exponential 
distribution of shower first interaction points, and then 
fit the entire distribution, and not just the trailing edge. 
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FIG. 1: B dependence on the pp total cross section app. 
The five curves are Unes of constant a''^^°!^-, of 414, 435, 456, 
499 and 542 mb — the central value is the published Fly's Eye 
value, and the others are ilcr and ±2(t. The solid curve is 
a plot of a QCD-inspired fit of B against app, obtained from 
a In^ fit — see text for details. The large dot is the prediction 
for at ^ = 77 TeV, the HiRes energy. 



as is done in the experiments of Ref. P, [3, 0, 13 • They 
obtain crP^f;^ = 460 ± 14 (stat) + 39 (syst) - 11 (syst) 
mb at y/s = 77 GeV, a result which they claim is ef- 
fectively model-independent and hence is an absolute 
determination fl3| . 

Extraction of app from crp'L°air- The total pp cross sec- 
tion is extracted from (7^™^^. in two distinct steps. First, 
one calculates the p-air total cross section, CTp_air, from 
the measured inelastic production cross section using 



prod 



q — c\ 



p— air 



(2) 



Next, the Glauber method |, 111] is used to transform the 
measured value of crp!°air i'^to a proton-proton total cross 
section app; all the necessary steps are calculable in the 
theory. In Eq. ([2]) the measured cross section for particle 



production is supplemented with CTp_air CpZair' 
elastic and quasi-elastic cross section, respectively, as cal- 
culated by the Glauber theory, to obtain the total cross 



section (Tp_air- 



The subsequent relation between a, 



prod 



and app critically involves the nuclear slope parameter 
B, the logarthmic slope of forward elastic pp scattering, 
daf /dt, i.e.. 



B 




(3) 



A plot of B against app, 5 curves of different values of 
is shown in Fig.[Tl taking into account inelastic 



prod 
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screening[l2j|. The reduction procedure from Cp^air to 



app is summarized in Ref. [8|. The solid curve in Fig. [T] is 
a plot of B vs. app, which we will discuss in detail later. 



Determination of app{s). Block and Halzen[l^ have 
made an analytic amplitude fit that saturates the Frois- 
sart bound p^. to both the available high energy total 
cross section and p-value data, where p is defined as the 
ratio of the real to the imaginary portion of the forward 
scattering amplitude, for both pp and pp interactions. 
For their high energy expressions they used the analytic 
amplitude form 



a^{v) = co + 



ciln(-) +C2ln2(-) +I3v' {-) 
\mJ \m/ Vm/ 



± 5 
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where the upper sign is for pp and the lower sign is for pp 
scattering, with /z = 0.5, v is the laboratory energy, /+(0) 
is a dispersion relation subtraction constant, and m the 
proton mass. The 7 real constants cq, ci, C2, /3-p' , 5, a and 
/+(0) are parameters of the fit. Since at high energies, 
s, the square of the cms energy, is given by 2mv, we see 
that their cross section approaches In^ s at high energies, 
thus saturating the Froissart bound[l3|. 

Using all of the cross sections, app and app, along with 
all of the p-values, ppp and ppp, in the Particle Data 
Group [15|] archive that were in the laboratory energy in- 
terval 18.3 <v< 1.73 X 10^ GeV, i.e., 6 < < 1800 
GeV, Block and Halzen[l^ formed a sieved data set us- 
ing the sieve algorithm of Ref. to eliminate outliers. 



which markedly improved their fit 13[ . Using 4 analyt- 
icity constraints 0, i.e., by fixing both the cross sections 
app and app and their laboratory energy derivatives, at 
~ 4 GeV, they reduced the number of parameters to 
be fit from 7 to 4 and obtained an excellent fit, which, in 
turn, constrained pp cross sections at cosmic ray energies 
to have a relative accuracy ~ 1 — 2%. Their fits to a and 
p are shown in Fig. [DJa) and Fig. [IJb), respectively. 

Determination of B(s). A QCD-inspired 

parameterization [3| of forward pp and pp scatter- 
ing amplitudes which is analytic, unitary and fits all 
data of CTtot, B and p for both pp and pp interactions 
has been made, using 2 analyticity constraints which fix 
app and app at ^/s = 4 GeV; see Fig.[2][c) for B{s). 

The solid curve in Fig. [1] is a plot of B vs. app, with 
B taken from the QCD-inspired fit of Re f. [ l6l| and app 
taken from the Froissart bound fit of Ref. [13| . The large 
dot corresponds to the value of app and B at ^/s = 77 
TeV, the HiRes energy, thus fixing the predicted value of 
^p-°air the HiRes energy. 

Obtaining app from o-p!°air- Fig. [31 we have plotted 
^prod ^Yiat are deduced from the 



the values of app vs. 



al 
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FIG. 2; The saturated Froissart bound fit[l3|| of total cross 
section app, p vs. ^/s, in GeV, for pp (squares) and pp (circles) 
accelerator data: (a) app, in mb, (b) p; (c) the nuclear slope 
B, in GeV~^ vs. ^/s, in GeV, from a QCD-inspired fitp^. 



intersections of the B-Upp curve with the crp™^^. curves in 
Fig.[T] Figure [3] furnishes cosmic ray experimenters with 
an easy method to convert their measured crpj_°^j. to Upp, 

and vice versa. The percentage error in crpl°^j. is ~ 0.4% 
near (7^^°^^. = 450mb, due to the error in app from model 
parameter uncertainties. 




FIG. 3: A plot of the predicted total pp cross section Opp, in 
mb vs. the measured p-air cross section, crpl°airi in nib. 

Determining the k value. It is important at this point 
to recall Eq. ([T]), A„j = fcAp_air, thus rewinding us of 
the fact that in Method I, the extraction of Ap_air (or 
CT^l°^j.) from the measurement of requires knowing 
the parameter k. The measured depth Xmax at which 
a shower reaches maximum development in the atmo- 
sphere, which is the basis of the cross section measure- 
ment in Ref. is a combined measure of the depth of 



the first interaction, which is determined by the inelastic 
cross section, and of the subsequent shower development, 
which has to be corrected for. The model dependent rate 
of shower development and its fluctuations are the origin 
of the deviation of k from unity in Eq. ([T|) . As seen in Ta- 
blclH its values range from 1.6 for a very old model where 
the inclusive cross section exhibited Feynman scaling, to 
1.15 for modern models with large scaling violations. 

Adopting the same strategy that earlier had been used 
by Block et al.[l3]j we decided to match the data to our 
prediction of o'pL°^j.(s) in order to extract a common value 
for k. This neglects the possibility of a weak energy de- 
pendence of k over the range measured, found to be very 
small in the simulations of Ref. 9]. By combining the 
results of Fig.[2](a) and Fig.[3l we obtain our prediction 
of cTp!"^!- vs. -^/s, which is shown in Fig. [D To deter- 
mine fc, we leave it as a free parameter and make a fit 
to rescaled <jI'°^.^{s) values of Fly's Eye, [i|AGASSA[3, 
EAS-TOP[| and YakutskQ, which are the experiments 
that need a common fc- value. 

Figure [4] is a plot of crpl°^j. vs. -^/i, the cms energy 
in GeV, for the two different types of experimental ex- 
traction, using Methods I and II described earlier. Plot- 
ted as published is the HiRes value at ^/s = 77 TeV, 
since it is an absolute measurement. We have rescaled 
in Fig.Hthe published values of 0-^"^^ for Fly's Eye[l|, 
AGASSA|2], Yakutsk[i] and EAS-TOPQ, against our 
prediction of (Tp™^j., using the common value oi k — 
1.264 ± 0.033 ± 0.013 obtained from a fit, and it is 
the rescaled values that are plotted in Fig. U) The er- 
ror in k of 0.033 is the statistical error of the fit, 
whereas the error of 0.013 is the systematic error due to 



the error in the prediction of cr. 
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p— air ■ 



Clearly, we have 



HiRes 

Yakutsk, rescaled 
FlysEye, rescaled 
AGASSA, rescaled 
EAS-TOP, rescaled 




FIG. 4: A X fit of the renormaUzed AGASA, EASTOP, Fly's 
Eye and Yakutsk data for 0-^1°^;^., in mb, as a function of the 
energy, ^/s, in GeV. The result of the fit for the parameter 
k in Eq. ((TJ is fc = 1.263 ± 0.033. The HiRes point (solid 
diamond), at -y/s = 77 GeV, is the model-independent HiRes 
experiment, which has not been renormalized. 
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an excellent fit, with complete agreement for all exper- 
imental points. Our analysis gave ^ = 3.19 for 11 de- 
grees of freedom (the low is likely due to overesti- 
mates of experimental errors). We note that our fc-value, 
k = 1.264 ± 0.033 ± 0.013, is about halfway between the 
values of CORSIKA-SIBYLL and CORSIKA-QSGSjet 
found in the Pryke simulations Q , as seen in Table [D 

We next compare our measured k parameter with a 
direct measurement of k by the HiRes group (3]. They 
measured the exponential slope of the tail of their 
distribution. Km and compared it to the p-air interac- 
tion length Ap_air that they found. Using Eq. ([1]), they 
deduced that k = 1.21 -I- 0.14 — 0.09, in agreement with 
our value, giving us additional experimental confirmation 
of our method. 

Conclusions. Our measured k value, k — 1.264 ± 
0.033 ± 0.013, agrees very well with the fc-value mea- 
sured by the HiRes group, at the several parts per mil 
level, and in turn, they both agree with Monte Carlo 
model simulations at the 5-10 part per mil level. 

It should be noted that the EASTOP0 cms energy, 
2 TeV, is essentially identical to the top energy of the 
Tevatron collider, where there is an experimental deter- 
mination of (Tpp[19|, and consequently, no necessity for 
an extrapolation of collider cross sections. Since their 
value of crpl°^j, is in excellent agreement with the pre- 
dicted value of CTp^^j,, this serves to anchor our fit at its 
low energy end. Correspondingly, at the high end of the 
cosmic ray spectrum, the absolute value of the HiRes ex- 
perimental value of tT^™^j. at 77 TeV — which requires no 
knowledge of the k parameter — is also in good agreement 
with our prediction, anchoring the fit at the high end. 
Thus, our CTp^jJj, predictions, which span the enormous 
energy range, 2 < < 100000 TeV, are completely 
consistent with all of the cosmic ray data, for both mag- 
nitude and energy dependence. 

In the future, we look forward to the possibility of con- 
firming our analysis with the high statistics of the Pierre 
Auger Collaboration [2^, as well as confirming the pre- 



diction of 107.3 ± 1.2 mb for the total pp cross section 13 1 
at the LHC energy of 14 TeV. 
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At the cms energy ^/s = 2 TeV, the total cross sections 
app and app are expected to differ by ~ 0.01 mb, and thus, 
for t he p urpose of predicting 0-^1°^^., are the same — see 
Ref. [g. 
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